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Abstract: 

We perform a Monte Carlo calculation of the second moment of the polarised valence 
quark distributions in the nucleon, using quenched Wilson fermions. The special feature 
of this moment is that it is directly accessible experimentally. At a scale of /i ~ 2 GeV 
we find = 0.198(8), A^d = —0.0477(33). We compare these numbers with recent 

experimental results of the SMC collaboration. 
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Measurements of polarised deep inelastic structure functions of the nucleon have re- 
vealed that only a small fraction of the nucleon's spin is carried by the spin of the 
quarks JI], 0, |3], f|, |^]. This raises the question as to how the spin of the nucleon is dis- 
tributed among spin and orbital angular momentum of its constituents. The constituents 
are valence quarks, sea quarks and gluons. In this paper we shall be concerned with the 
spin carried by the valence quarks. 

In the naive quark parton model the unpolarised structure functions and the polarised 
structure function gi can be expressed in terms of the probability distributions for finding 
a quark with spin parallel, q^, and antiparallel, qi, to the longitudinally polarised nucleon: 

Fi{x) = Wel[q(x)+q(x)\, (1) 
z q 

9i(x) = We 2 g [Aq(x) + Aq(x)), (2) 

Z q 

where 

q(x) = g T (x) + q l (x), Aq(x) = g T (x) - q^x). (3) 
By analogy one can define hadron inclusive structure functions 

F^z) = \Y,^) D >) + ^) Dh q^)l ( 4 ) 
/ q 

fit?**) = ^%fa{x)Di(z)+Aq(x)Di(z)), (5) 

Z q 

where D^(z) is the fragmentation function for quark q to produce hadron h, and z = 
p h -p/p-q with p h , p and q being the hadron, nucleon and photon momentum, respectively. 
Of interest to us is the polarisation asymmetry M 

A h _ 9ii^, z ) / fi x 

Consider the polarisation asymmetry of tt + minus tt~ inclusive cross sections. For a proton 
and deuteron target we find 

A7V+ _ W - = 4Au™\x) - Ad™\x) 
v 4u val (x) - d val (x) 

and 

+ __ Au^(x) + Ad™ l (x) 

d u™\x) + d™\x) ' 1 ' 

respectively. Note that the fragmentation functions drop out completely and that the sea 
quarks do not contribute. This happens because of isospin invariance relating the various 
fragmentation functions with each other. The advantage of these quantities is that they 
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allow a direct measurement of the valence quark distribution functions || . The latter can 
then be compared with the results of quenched lattice calculations J7|. 

The polarisation asymmetries (|7|), (§) have been measured by the SMC collaboration. 
For the first moment 

Aq = / dxAq(x), q = u,d (9) 



o 

it has been found, [|, Aw val = 1.01 ± 0.19 ± 0.14 and Ad val = -0.57 ± 0.22 ± 0.11. The 
lattice results are, 0, Au = 0.841(52) and Ad = -0.245(15), which include data on 24 3 32 
lattices. The analysis has recently been extended to the second moment 

A (1) g = [ dxxAq(x). (10) 



o 



The SMC collaboration has obtained results for A^u and A^d [T(|. A similar analysis 
will be published soon by the HERMES collaboration, [±T|. 

In QCD, A«g = \a!f where a? is the matrix element of the operator 



0^=( t -)^^D iM> q t (11) 



defined by 



(p, s\SO { 5 M) ^\p, s) = a[ g) Ss {M ^p iM)v , (12) 



where S symmetrises and subtracts traces. We have emphasised Minkowski space with 
an index Ai. s^ M ' is the nucleon spin vector defined by 

» ( M¥4 '=™{ g +^hf)' (i3) 

with e being the rest frame quantisation axis and o = ±1 the spin component along this 
axis. We now euclideanisef] and discretise the operator given in eq. (|TI|). Let us first 
define our euclidean operator by 

O^u = Du q. (14) 

Its transcription to the lattice is straightforward. However on the lattice the symmetry 
group is reduced from 0(4) to the hypercubic group H(A). This loss of symmetry increases 
the possibilities of mixing under renormalisation. In the continuum, application of S is 
sufficient to construct operators which are multiplicatively renormalisable (in the flavour- 
nonsinglet sector). To achieve the same on the lattice one has to work slightly harder 



T3H . In our case one finds two multiplets of operators (corresponding to two inequivalent 



irreducible representations of H{A)) which are multiplicatively renormalisable. They are 
1 Conventions are given in ref. 
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0(a) 


(p, s\O u \p, s) 


Representation 




1 < H < v < 4 
1 < i < 3 


-aiSiPi 


rj 6 \C = -l 
rf,C = -l 



Table 1: Multiplicatively renormalisable lattice operators 0(a) and the continuum matrix 
elements of the corresponding renormalised operators. {• • •} means symmetrisation. 
denotes the representation, 



13], with m being its dimension. C is the charge conjugation 



eigenvalue of the operator. 



given in table [1]. For each of them we obtain the renormalised operator R (fi) (renormal- 
isation scale fi) from the lattice operator 0(a) (lattice constant a) by multiplying with 
the appropriate renormalisation factor Zq\ 

R (ii) = Z (M 2 ,g(a))0(a). (15) 

In our computation we have only considered p = or p = p\ = p\e\ together with 
spin-quantisation axis e = e 2 ■ Thus from Table [I] we can consider 

O ai>a = 05{12}, 

O au b = 5{42 }- (16) 

Note that for a non-zero matrix element for O ai jQ we must have a non-zero three-moment- 
um, while for O ai ^ we do not have this restriction. As including three-momentum in the 
matrix element makes the signal more noisy we would expect the best result to come from 
using zero three- momentum together with the matrix element of O ai ,b- 

The nucleon matrix elements are computed from ratios of three- to two- point func- 
tions. Thus with 

C r (t,p) = Y, r ^(B a (t,p)Bp(0,p)), 

a/3 

C r (t,T-p,O ai ) = J2 r 0c(M-t,P)O ai (r)Bp(O,p^), (17) 

a/3 

and 

R(t,r;p,O ai ) 

we find 

R(t,T;p,O aua ) 
R(t,r;p,O aub ) 



Ci n , .. (t, r;p, O ai ) 
2(1+74)47572 v ^ 1J 



1 1 ttin 

P\ai t > r > 0, 



Z Qai 2kAE n 



Z 2k 4 



1 % -m N a x t > r > 0, (19) 



3 



where k is the Wilson hopping parameter and tjin {En) is the mass (energy) of the 
nucleon. Eq. [19| holds provided the operator is inserted at a Euclidean time r sufficiently 
far away from the nucleon source at time t — and the sink at time t. We then expect 
to observe a 'plateau' where R is independent of r. 

We have simulated quenched Wilson fermions with five different k values on iV s 3 x 32 
lattices with N s = 16, N s = 24 at (3 = Q/g 2 = 6.0. p\ is taken to be the lowest possible 
momentum, namely p\ = 2tt/N s . More details of the computation may be found in 
ref. 0. Note, in particular, that we have dropped the quark-line disconnected terms. 
This, coupled with the use of the quenched approximation, means that sea quarks play 
only a small role in our calculation, so that we are effectively measuring valence moments 
only. The masses for this calculation are taken from |15| . In Fig. |l| we show typical ratios 
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Figure 1: The R = R(t = 13, r,p, O ai ) ratio at k = 0.1530. 

for t = 13. The fit interval for the determination of the plateau is chosen to lie between 
r = 4 and r = 9. These values of r seem to be sufficiently far away from both the baryon 
source t = and sink t = 13 to avoid contamination with excited states. 
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To find A^u and A^d, we must also calculate the renormalisation constant Z^f . 
This was given in ref. |T2| where it was called Zq t , see also [pj] . Explicitly, we have 



Z^(M\ g) = l- ^_C F [ lQai Info*) + + 0(g 4 ), (20) 

with Cf = 1,70a! = y, an d -Sol 5 — 0.3454. We shall work, for convenience, at a 
scale of n = 1/a. This gives for (3 = 6.0, Z 0oi (l,l) = 0.997. We see that the 0(# 2 ) 
coefficient is thus very small. This is because 'tadpole terms' in the perturbation expansion 
have cancelled. Indeed if we try to improve the result using 'tadpole improved boosted 
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perturbation theory', [18], we find 



^(M 2 , <?)) = ^ (l - ^^ C ^o ai InM + B*£ + (n D - l)^ 2 ] ] , (21) 



i 



where rip is the number of derivatives, Uq = (^tr[/ n )4 with Un the expectation value of 
the plaquette, and a MS (ii) is used here as the boosted coupling. In our case n D = 1, so 
the coefficient of the coupling constant is unchanged and remains small, and at /3 — 6.0, 
a m (l/a) w 0.1981, 0, so that 



1) = 0.993 (22) 

which is a negligible change in the previous result. Our conclusion is that the uncer- 
tainty in the renormalisation constants is probably quite small, although only a full non- 
perturbative calculation would reveal this of course. We shall use, in the following, the 
result from eq. (|22|) . 

Finally we note that to determine the scale [l for our results we must first estimate 
a. Measuring a physical quantity on the lattice (such as a mass) and comparing this to 
the experimental value determines a. However as Wilson fermions have 0(a) corrections, 
it is better to set the scale from a gluonic quantity, as this has only 0(a 2 ) corrections. 
This also avoids the need of making a chiral extrapolation in the masses. From the string 
tension K we get (using a scale of \/K = 427MeV, [0)0 

a' 1 = 1.95(2) GeV, (23) 

where the error is purely statistical. 

We now come to our results. In Tables |2], [| we give our numbers and in Fig. |2] we 
plot A^m, A^d against the bare quark mass 

om « = \ (; " z) • (24) 



2 From a combined fit of lattice data, see ref. |fl5| . 
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K 


0.1515 


0.1530 


0.1550 


# configs 


0(400) 


0(500) 


0(900) 


u \P=Pi 

A (1) ?;l - 

u \p=Pi 
u \p=Pi 

^6 u lp=pi 


23fif1 3"l 
0.241(10) 
0.244(18) 
-0.0577(52) 
-0.0596(27) 
-0.0596(48) 


940(1 X) 
0.234(9) 
0.229(19) 
-0.0577(60) 
-0.0559(24) 
-0.0586(51) 


21 Qf2fT) 
0.223(10) 
0.233(24) 
-0.0549(148) 
-0.0556(38) 
-0.0651(10) 



Table 2: The results for A^q on a 16 3 x 32 lattice. 



K 


0.1550 


0.1558 


0.1563 


# configs 


O(100) 


O(100) 


O(100) 


u \p=pi 

A W u\ - 

^b a \p=0 
^6 u \p=pi 
u \p=Pi 

A [1) d\^ - 

^6 u lp=Pi 


0.200(25) 
0.219(14) 
0.223(21) 
-0.0437(81) 
-0.0526(48) 
-0.0559(63) 


0.172(32) 
0.211(14) 

0.212(25) 
-0.0409(141) 
-0.0463(68) 
-0.0473(95) 


0.122(58) 
0.193(17) 
0.184(37) 
-0.0755(293) 
-0.0416(103) 
-0.0344(195) 



Table 3: The results for A^q on a 24 3 x 32 lattice. 

together with a linear chiral extrapolation. (The chiral limit is determined from the 
vanishing of the pion mass, which occurs here at k c = 0.15721(1), ||T5| . We roughly 
estimate our quark masses from am^ = Z m am q ~ 1.21ara 9 , fllS] , to be m MS (2GeV) ~ 
290, 210, 110, 70, 45 MeV.) The quality of the data seems quite good. The statistics 
for the light quark masses is rather low, so although we include them in the linear fit, 
they should be regarded as only confirming the heavier quark mass results. Note that the 
results for k = 0.1550 have been measured on two different lattice sizes and differ by less 
than their errors which would indicate that at least up to this value of k finite size effects 
are under control. At each k we have three different measurements of A^q, one at p = 
and two aXp — p\. If Lorentz invariance has been restored, all three measurements should 
agree, which is what they do within errors. As expected the p = results are the best, 
although the other signals are quite acceptable. It would also seem that the gradient in 
the chiral extrapolation is rather small over this quark mass range. The values in the 
chiral limit are given in Table |j. 

For definiteness, averaging the three different measurements we quote a result of 

A (1) m = 0.198(8), 

A {1) d = -0.0477(33), (25) 
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Figure 2: The chiral extrapolation for A^q. 

at a scale of /i 2 m 4 GeV 2 . 

We first compare our result with that given in pOf where a phenomenological fit was 



made from the experimental data for the polarised structure function g± to obtain the 
polarised parton distribution functions. Taking the second moment of their NLO fit 
functions yields values of around A^-'u ~ 0.15, A^d m —0.055, at a scale of /x 2 = 4 
GeV 2 . While the d moment is in agreement with our result, their u moment is a little 
smaller than ours. 

We can see how the matrix element runs with the scale upon using the renormalisation 
group equation 



c F 



A w M = l^±\ 2b0 A (l) ( ) (26) 

with the value of a m (2GeV) « 0.1981, (T§, and with b = 11 - |n/ = 11, as we are 
working in the quenched approximation. Thus at /i 2 = lOGeV 2 we find: 

A (1) m = 0.189 ±0.008, 
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A^u\- - 

a \P=Pi 




A (1) ?/l- - 

^b a \P=Pi 






Afe Mp=Pi 


0.176(23) 


0.201(10) 


0.203(20) 


-0.0427(94) 


-0.0477(41) 


-0.0507(71) 



Table 4: The results for A^ q in the chiral limit. 

A {1) d = -0.0455 ±0.0032. (27) 

No great dependence on the scale is seen. 

Finally our predictions can now be compared with the preliminary results of the SMC 
collaboration ]10| at the scale fi 2 = 10 GeV 2 which read: 



A w u = 0.169 ±0.018 ±0.012, 

A w d = -0.055 ±0.027 ±0.011, (28) 

where the first quoted error is statistical and the second one systematic. The integrals have 
been evaluated using data from the measured range 0.003 < x < 0.7. The unmeasured 
regions below x = 0.003 and above x = 0.7 give probably a negligible contribution to A^u 
and A^d. The lattice and experimental results agree within their respective errors. 



ACKNOWLEDGEMENTS 

The numerical calculations were performed on the APE (Quadrics QH2) at DESY-Zeuth- 
en, with some of the earlier computations on the Bielefeld University APE. We thank 
both institutions for their support. One of us (LM) was supported by BMBF and the 
KBN grant 2 P03B 065 10. MG and PR gratefully acknowledge support by the Deutsche 
Forschungsgemeinschaft . 



References 

[1] M.J. Alguard et al., Phys. Rev. Lett. 37, 1258 (1976); ibid 37, 1261 (1976); 
G. Baum et al., Phys. Rev. Lett. 51, 1135 (1983). 

[2] J. Ashman et al., Phys. Lett. B206, 364 (1988); 
J. Ashman et al., Nucl. Phys. B238, 1 (1989). 

[3] B. Adeva et al., Phys. Lett. B302, 500 (1993); ibid B320, 400 (1994); 
D. Adams et al., Phys. Lett. B329, 399 (1994); ibid B357, 248 (1995). 

[4] P. L. Anthony et al., Phys. Rev. Lett. 71, 959 (1993). 

[5] K. Abe et al., Phys. Rev. Lett. 74, 346 (1995); ibid 75, 25 (1995); 
K. Abe et al., Phys. Lett. B364, 61 (1995). 



8 



[6] L. L. Frankfurt, M. I. Strikman, L. Mankiewicz, A. Schafer, E. Rondio, A. Sandacz 
and V. Papavassiliou, Phys. Lett. B230, 141 (1989). 

[7] M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz and A. 
Schiller, Phys. Rev. D53, 2317 (1996). 

[8] B. Adeva et al., Phys. Lett. B369, 93 (1996). 

[9] C. Best, M. Gockeler, R. Horsley, L. Mankiewicz, H. Perlt, P. Rakow, A. Schafer, 
G. Schierholz and S. Schramm, preprint DESY 97-116 (1997) flhep-ph/9706503 ), to 
appear in the proceedings of DIS97. 

[10] J. Pretz, Dissertation, Mainz (1997). 

[11] M. Diiren, private communication. 

[12] C. Best, M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, A. Schafer, 
G. Schierholz and S. Schramm, preprint DESY 97-41 (1997) ( |tiep-lat/9703014| ), to 
appear in Phys. Rev. D. 

[13] M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz and A. 
Schiller, Phys. Rev. D54, 5705 (1996). 

[14] M. Gockeler, R. Horsley, H. Perlt, P. Rakow, G. Schierholz, A. Schiller and P. 
Stephenson, Phys. Lett. B391, 388 (1997). 

[15] M. Gockeler, R. Horsley, H. Perlt, P. Rakow, G. Schierholz, A. Schiller and P. 
Stephenson, preprint DESY 97-125, HUB-EP-97/38 (1997) ( |hep-lat/9707021] ). 

[16] M. Gockeler, R. Horsley, E.-M. Ilgenfritz, H. Perlt, P. Rakow, G. Schierholz and A. 
Schiller, Nucl. Phys. B472, 309 (1996). 

[17] M. Gockeler et al., in preparation. 

[18] G. P. Lepage and P. B. Mackenzie, Phys. Rev. D48, 2250 (1993). 

[19] E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane and T. M. Yan, Phys. Rev. D21, 
203 (1980). 

[20] T. Gehrmann and W. J. Stirling, Phys. Rev. D53, 6100 (1996). 



9 



